Quantitative prediction of the phase diagram of DNA-functionalized nano-colloids 
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We present a coarse-grained model of DNA-functionalized colloids that is computationally 
tractable. Importantly, the model parameters are solely based on experimental data. Using this 
highly simplified model, we can predict the phase behavior of DNA-functionalized nano-colloids 
without assuming pairwise additivity of the inter-colloidal interactions. Our simulations show that 
for nano-colloids, the assumption of pairwise additivity leads to substantial errors in the estimate 
of the free energy of the crystal phase. We compare our results with available experimental data 
and find that the simulations predict the correct structure of the solid phase and yield a very good 
estimate of the melting temperature. Current experimental estimates for the contour length and 
persistence length of single-stranded DNA sequences are subject to relatively large uncertainties. 
Using the best available estimates, we obtain predictions for the crystal lattice constants that are 
off by a few percent: this indicates that more accurate experimental data on ssDNA are needed to 
exploit the full power of our coarse-grained approach. 

PACS numbers: 07.05.Tp, 64.70.Nd, 64.75.Yz, 82.70.Dd 



Nature provides spectacular examples of complex, 
functional systems that self-assemble from small, pre- 
fabricated units. The grand challenge in nano-material 
design is to imitate this phenomenon to construct com- 
plex, functional materials. To achieve this goal we must 
be able to design suitable nano-sized building blocks 
and to create conditions that cause these entities to 
self- assemble into the desired target structures. DNA- 
based building blocks — which take advantage of the se- 
lective binding of bases on complementary strands to 
guide assembly — offer attractive model systems to ex- 
plore self-assembly strategies (see e.g. Refs. [E Hi)- For 
instance, colloidal particles can be functionalized with 
a short single-stranded (ss) DNA sequence tethered to 
an inert polymeric "spacer" 0]. These "sticky ends" on 
two different colloids may then either bind directly to 
each other via complementary sequences or via a ssDNA 
linker sequence introduced in solution, allowing the col- 
loids to form three-dimensional structures. The proper- 
ties of these building blocks, temperature, pH and ionic 
strength of the parent solution determine if, and on what 
time scale, self-assembly takes place. 
In view of the vastness of this "design space" [4j, a 
careful selection of experimental conditions is crucial: 
experiments on nano- and micron-sized DNA- 



functionalized colloids (DNACs) [llj, |12[ show that self- 
assembly of spatially ordered structures requires consid- 
erable fine-tuning: under most experimental conditions, 
amorphous aggregates form, even if the ordered (crys- 
tal) phase is thermodynamically stable. The experiments 
show that the stability of crystals depends crucially on 
temperature (T) and on the length and flexibility of the 



spacers. Of the crystal structures experimentally ob- 
served, the bcc and fee lattices are the most stable, while 
to date little progress has been made with self-assembly 
of DNAC crystals that are more complex than cubic. 
Therefore, it is desirable to have theoretical and nu- 
merical guidance in selecting optimal conditions for self- 
assembly. However, an "ab-initio" approach to study the 
phase behavior of DNACs quantitatively, using many- 
particle simulations, goes well beyond the state-of-the- 
art. Existing theoretical studies are based on the as- 
sumption that the interaction between DNACs is pair- 
wise additive [l3[ , and numerical simulations of DNACs 
typically employ highly simplified, ad-hoc coarse-grained 
models that also assume pairwise additivity of interac- 
tions (see e.g. Refs. l4-|22[). An exception is the work 
of Rogers et al. [H, QA\. However, the approach followed 
in their paper, whilst suited to compute pair potentials 
between micro-sized DNACs, would become computa- 
tionally prohibitively expensive for nano-sized DNACs 
where, as we will show below, it is not warranted to 
assume pairwise additivity of the inter-colloidal interac- 
tions. What is more, those systems that have been shown 
most promising for crystallization, are nano-DNACs cov- 
ered with DNA strands in an intermediate length regime, 
where the DNA strands are too long to be approximated 
within the rod-picture [20|, but too short to make use of 
polymer scaling laws. As a consequence, existing models 
cannot be used to predict the stability of crystal phases 
of nano-DNACs. 

In this Letter, we present a quantitative numerical ap- 
proach to predict the thermodynamic stability and phase 
behavior of nano-DNACs based on a coarse-graining pro- 
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cedure free of fitting parameters. 

We validate our coarse-graining approach on a system 
that has been studied extensively in experiments @, 
namely a symmetric, binary mixture of gold nano-colloids 
of a radius Rc ~ 6nm which were grafted with ~ 60 ss- 
DNA strands. Colloidal species A and B only differ in the 
sequence of their sticky ends, which mediate the binding 
between colloids A and B. In the experiments of Ref. [6j, 
it was observed that systems with ssDNA strands of more 
than 50 nucleotides crystallize into CsCl structures. In 
our study, we therefore focus on their system of DNACs 
functionalized with ssDNA strands of 65 nucleotides, 15 
of which were responsible for binding. These DNACs 
were experimentally found to crystallize for temperatures 
below 7^=62. 5°C. 

In order to arrive at a computationally tractable model, 
we carry out a staged coarse-graining procedure. At the 
most microscopic (yet not fully atomistic) level, we rep- 
resent the Au colloid as hard sphere of radius Rc and 
we model the ssDNA strands as freely jointed, charged 
chains with a Kuhn length of 1.5nm 25j, [26J, and using 
an inter-base distance for ssDNA of 0.43nm 26]. We 
stress, however, that the values reported in these (and 
other) papers are subject to considerable error bars, and 
are likely to be sensitive to both the physical conditions 
of the solution and the precise DNA sequence. Hence, 
we should expect that these inaccuracies will translate 
into errors in computed characteristic length scales, such 
as the lattice constants of DNAC crystals. As we show 
below, this is precisely what we observe. 
The charge carried by the ssDNA's backbone is assumed 
to be distributed equally among the vertices of the freely 
jointed chain. The vertices interact through a Debye- 
Hiickel potential that depends on T, and on the dielec- 
tric constant and ionic strength of the solvent j25|. We 
assume that the grafted DNA is distributed uniformly on 
the Au surface and does not diffuse [27j |. 
In spite of the simplicity of this model, it has a large num- 
ber of degrees of freedom per colloid (~ 3000). As a con- 
sequence, simulations are only feasible for relatively small 
systems. In order to be able to treat DNA hybridiza- 
tion in systems containing many colloids, we therefore 
develop the next level of coarse-graining: a "core-blob" 
model in which each sticky end is coarse-grained to a 
polymer "blob" , while the bare colloid and the remaining 
segments of all strands constitute an effective colloidal 
"core" . In this way, we reduce the degrees of freedom 
per colloid tenfold. We determine the parameters and in- 
teraction potentials that characterize the core-blob model 
in Monte Carlo simulations of the underlying microscopic 
model [28] . Further, we allow for DNA hybridization in 
the core-blob model. The binding probability of comple- 
mentary DNA strands is computed on the basis of the 
tabulated hybridization free energy of two complemen- 
tary ssDNAs in solution (29| . which is directly based on 
experimental data. Details of the coarse-graining pro- 
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FIG. 1. (Color online) Left: The pair interaction $£ B (r,T) 
(bold lines) is the sum of the T-independent steric repulsion 
&2,re P (r) (black thin line) and the strongly T-dependent at- 
tractive hybridization energy 3?2,hyb(r, T) (colored thin lines). 
Right: The three-body effective interaction /3$3" BA (r) (bold 
solid lines) of two A and one B colloids arranged in an equi- 
lateral triangle (see inset) compared with the sum of the two- 
body contributions 3/3$2, rC p + 2/3$2,hyb (bold dashed lines). 
Further, the repulsive and the attractive hybridization con- 
tributions to the interactions are shown explicitly in thin 
lines (two-body: dashed lines, three-body: solid lines). Both: 
Results are shown at 56.9°C, 59.1°C, 61.4°C, 63.2°C, and 
65.1°C. Inset: A simulation snapshot of 1 B colloid (red solid 
sphere) interacting with 2 A (green solid spheres) at 59.1°C 
and r — 5.75Rc- Free blobs are drawn as yellow spheres, 
while hybridized sticky ends are shown as blue rods. The 
translucent spheres indicate the extension of the core. 



cedure and the Monte Carlo scheme used to carry hy- 
bridization moves will be presented elsewhere [30| . 
In what follows, we will use the core-blob model to pre- 
dict the phase behavior of DNACs. However, we first 
go one step further in coarse-graining and compute the 
effective pair potentials ^(r, T) between two DNACs 
(H = {AA,BB,AB}) as a function of the distance r and of 
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(r, T) accounts for the steric repulsion $ 2 



rep 



T. $ 

between the colloids, which can be determined as de- 
tailed in 28[. In the AB case, there is an additional 
attractive hybridisation interaction < I > 2,hybj which is pro- 
portional to the amount of hybridized DNA strands and 
is calculated by thermodynamic integration as described 
in Ref. [2(j. We find that $2. rep is nearly constant be- 
tween 25 and 75°C, while $2,hyb is strongly T-dependent 
(Fig.Q]). For T < X^ xp , $ AB therefore develops a strongly 
T-dependent minimum at distances r ~ 5.25 to 6i?c- For 
the system presented here, a T-difference of only 5.2°C 
(62.1 to 56.9°C) leads to a difference in minimum in $ AB 
of ~ 20ksT ~ 13kcal/mol. The difficulty to make high- 
quality crystals of DNACs is related to this strong tem- 
perature dependence of $ 2 VB : defects in growing crystals 
can only anneal if T is just below the melting tempera- 
ture. At lower T, bonds are very stable and dense aggre- 
gates, once formed, cannot equilibrate. 
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We are unaware of direct measurement of the pair po- 
tentials of DNA-functionalized nemo-colloids, hence we 
cannot validate our model at this level. However, nu- 
merically, we can test whether the assumption of pair- 
wise additivity of $2 > as typically employed in theoretical 
studies, is warranted. To this end, we compute the inter- 
action $|f(r,T) (N = {AAA, BBB, ABA, BAB}) between 
a triplet of colloids arranged on the vertices of an equi- 
lateral triangle with sidelength r (inset Fig. [T]) . We find 
that <j> AAA > BBB — $ 3 , rcp is well represented by 3$2,rc P for 
all T considered. For the ABA (=BAB) case, additivity 
would imply: <t>$ BA {r,T) = 3$ 2 , rep (r,T) + 2$ 2)hyb (r,T). 
As can be seen in the right panel of Fig. [TJ even for tem- 
peratures close to T™ p this relation is only fulfilled at 
large separations r > 6Rc- The reason why additivity 
fails in this case is that different B colloids compete for 
the same DNA on A. This effect is not accounted for in 
$2, which measures the average amount of bonds in an 
isolated AB pair and hence overestimates the attraction 
between AB pairs in a trimer 31 1. 

Turning our attention to multi-particle systems, we 
study the thermodynamic stability of various crystalline 
phases. Using the pair potentials, we first identify can- 
didate equilibrium crystal structures via a genetic algo- 
rithm similar to the one used in Refs. [32I |33|. This ap- 
proach, which neglects positional entropy, suggests that 
the CsCl structure is most stable close to the melting 
density, while a NaTl structure is predicted to be sta- 
ble at higher densities. Experimentally, the distinction 
between these two structures is not straightforward as 
X-ray scattering only probes the arrangement of the Au 
cores which is the same for CsCl and NaTl crystals. We 
also considered CuAu, NaCl, 'straight' hep (s-hep), ZnS, 
AuCd and substitutional^ disordered CsCl and CuAu 
structures. For all crystal structures and also the fluid 
phase, we compute the free energy F for T < T™ p via 
thermodynamic integration within the core-blob model 
(Fig. [2]). The CsCl phase is found to have the lowest 
free energy for low colloidal volume fractions r). Only 
slightly higher in F, we find metastable s-hep, CuAu and 
NaTl structures. While the CsCl is already mechanically 
stable for 77 = 0.023, the latter structures are only me- 
chanically stable for slightly higher 77 > 0.026 — 0.035. 
At high 77 ~ 0.065 a metastable AuCd phase appears. 
The NaCl and ZnS structures are mechanically unstable 
for all r\ considered. For 77 > 0.07, the NaTl structure 
(composed of two interpenetrating diamond structures) 
competes — as predicted — with the CsCl structure; how- 
ever, this is in a density regime where the validity of our 
model is not guaranteed. Moreover, such dense crystals 
cannot easily be prepared in experiments. 
We stress that the formation of the low-density crystals is 
due to the DNA links between colloids. At much higher 77 
the crystal structure will be dictated by excluded volume 
interactions rather than DNA links. Then, close-packed 
structures such as CuAu and s-hep should be more stable 
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FIG. 2. (Color online) The dimensionless free energy PnF/N 
at T = 56.9° C as function of the colloidal volume fraction 77 
shows that the CsCl structure is the most stable one. Dashed 
lines: pair potential approach; solid lines: core-blob approach. 
CsCl: A, s-hcp:«, CuAu:B, NaTl:*, disordered CsCl:*, disor- 
dered CuAu:x, liquid: -4. The common tangents between the 
equilibrium CsCl crystal and the dilute vapor are shown as 
dotted lines for both models. Right inset: the lattice con- 
stant a as function of T in °C as measured in the experiments 
(heating: A, cooling: T; data from Ref. @]) and as obtained 
from simulations with the core-blob model (•) and using the 
pair potential approach (■). The data of a have been scaled 
to the values of a Ic f — a(56.9°C) as measured with the re- 
spective approaches: a°* f p = 45.3nm, a corc-biob _ 39 7 nm anc j 

a*j = 36.9nm. In the shaded region crystals are not stable 
in experiments. Left inset: a cut through a simulated CsCl 
structure at 77 = 0.024 and T = 56.9° C. For better visualiza- 
tion, only hybridized links are drawn. 



than the more open structures that dominate at lower 77. 
We also consider substitutional disorder and find that 
strongly disordered crystal structures are only mechan- 
ically stable for 77 > 0.05 and less stable than the cor- 
responding ordered structures. However, some substitu- 
tional disorder is inevitable: at 56.9°C, the free-energy 
cost of a single AB exchange in an otherwise ordered 
crystal structure ranges from AF = O.bksT (77 = 0.024) 
to 1M B T (77 = 0.05) for CsCl and from AF = 0.1fc B T 
(77 = 0.026) to 1M B T (?/ = 0.05) for CuAu. We note 
that substitutionally disordered crystals have been ob- 
served both in experiment and in simulations [22j ]. 
The present results suggest that slightly disordered struc- 
tures are metastable but kinetically arrested for T < T™ p 
while strongly disordered crystals will melt at low 77. 
In Fig. [21 we further compare F as obtained from the 
core-blob approach with the corresponding quantity ob- 
tained employing the pair potentials. The assumption of 
pairwisc additivity of $2 leads to a serious under-estimate 
of F. Nevertheless, for T < T°f p this approach provides 
a reasonable estimate of the range of mechanical stabil- 
ity of the various crystal structures and predicts the same 
phase order as the core-blob model, albeit in a more com- 
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pressed r\ range. 

We also used the core-blob model to test whether the 
colloidal suspension could undergo a transition between 
a dilute and a concentrated unordered phase below the 
freezing density and for T < T^ xp . We find that this is 
not the case: the CsCl-structure directly coexists with a 
very dilute suspension [34|. Using the common-tangent 
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construction (see Fig. [2J , we can determine the volume 
fraction — and thereby the lattice constant a — of the CsCl 
structure at coexistence. Comparison of the experimen- 
tally determined a at various T with the values obtained 
from the core-blob approach (right inset Fig. [2]) shows 
that the latter approach can predict the correct tem- 
perature dependence of the lattice spacing. This is not 
the case for the pair potential approach, which predicts 
an excessive contraction of the crystals. Both compu- 
tational approaches predict lower values of a than ob- 
served experimentally. The core-blob approach under- 
estimates a by ~12%, a discrepancy to be expected in 
view of the incomplete experimental information on the 
contour length and persistence length of ssDNA. With 
better experimental data, we expect that the core-blob 
model will yield a quantitative prediction of the DNAC 
crystal lattice constants. The predictions of the lattice 
constants as obtained by use of pair potentials are qual- 
itatively wrong - in particular this approach incorrectly 
predicts the temperature dependence of the lattice spac- 
ing. We can further compare the experimentally deter- 
mined crystal melting point with the melting point of 
simulated crystal slabs in equilibrium with the dilute va- 
por. In these simulations, we find that crystals melt for 
T > 64.3(5)°C, which matches well with the experimen- 
tally determined melting temperature X^ xp = 62.5(3)°C 
for the 65 nucleotide system of Ref . @ . The good quan- 
titative agreement between simulation and experiment in 
melting temperature is significant as the core-blob model 
contains no a posteriori adjustable parameters. 
The present approach has been developed to model 
DNACs with direct hybridization between sticky ends. 
However, our approach can be generalized to describe 
binding via linker sequences and it can be adapted to 
describe asymmetric and polydisperse systems, to name 
but a few examples. Although some adjustment of the 
underlying microscopic model of DNACs is needed to 
mend the accuracy of length predictions, our approach 
offers a path to the computer-aided design of DNA- 
functionalized building blocks that could be used to con- 
struct truly complex self-assembling structures. 

We thank B. Capone, S. Angiolctti-Uberti, 
B. M. Mognetti, W. Jacobs and G. Day for useful 
discussions. BMM acknowledges EU funding (FP7- 
PEOPLE-IEF-2008 No. 236663) and the MFPL VIPS 
program (funded by the BMWF and City of Vienna). 
AD was supported by an EMBO longterm fellowship. 
DF and FJMV acknowledge support of ERC Advanced 
Grant 227758. DF acknowledges a Wolfson Merit Award 



[1 
[2 
[3] 
[4] 
[5 
[6 
[7 

[8 
[9 

[10 
[11 

[12 
[13 
[14 

[15; 
[16 

[17; 

[is; 

[19 

[20 

[21 
[22 
[23 

[24 

[25 
[26 
[27 
[28 
[29 

[30 

[31 



[32] 
[33] 
[34] 



bianca.mladek@univie.ac. at 
df246@ca m.ac.ukl 

R. Bashir, Supperlattice Microstr. 29, 1 (2001). 
F. A. Aldaye et al, Science 321, 1795 (2008). 

C. A. Mirkin et al, Nature 382, 607 (1996). 

N. Geerts and E. Eiser, Soft Matter 6, 4647 (2010). 
S. Y. Park et al, Nature 451, 553 (2008). 

D. Nykypanchuk et al, Nature 451, 549 (2008). 

R. J. Macfarlane et al, Proc. Natl. Acad. Sci. 106, 10493 

(2009) . 

H. Xiong et al, Phys. Rev. Lett. 102, 015504 (2009). 
R. Macfarlane et al, Angew. Chem. Int. Ed. 49, 4589 

(2010) . 

P. Cigler, et al. Nat. Mater. 9, 918 (2010). 

P. L. Biancaniello et al, Phys. Rev. Lett. 94, 058302 

(2005). 

A. J. Kim, et al. Langmuir 22, 1991 (2006). 

A. V. Tkachenko, Phys. Rev. Lett. 89, 148303 (2002). 
D. B. Lukatsky and D. Frenkel, Phys. Rev. Lett. 92, 
068302 (2004).' 

F. Pierce et al, Langmuir 21, 8992 (2005). 

F. W. Starr and F. Sciortino, J. Phys.: Condens. Mat. 

18, L347 (2006). 

B. Bozorgui and D. Frenkel, Phys. Rev. Lett. 101, 045701 
(2008). 

R. Dreyfus, et al. Phys. Rev. Lett. 102, 048301 (2009). 
F. J. Martinez- Veracoechea, et al. Soft Matter 6, 6136 
(2010). 

M. E. Leunissen and D. Frenkel, J. Chem. Phys. 134, 
084702 (2011). 

R. T. Scarlett et al, Soft Matter 7, 1912 (2011). 

C. Knorowski et al, Phys. Rev. Lett. 106, 215501 (2011). 
W. B. Rogers and J. C. Crocker, Proc. Natl. Acad. Sci. 
108, 15687 (2011). 

B. M. Mognetti et al, Proc. Natl. Acad. Sci. 109, E378 
(2012). 

Y. Zhang et al, Biophys. J. 81, 1133 (2001). 

B. Tinland et al, Macromolecules 30, 5763 (1997). 

S. Stranick et al, J. Phys. Chem. 98, 11136 (1994). 

B. M. Mladek and D. Frenkel, Soft Matter 7, 1450 (2011). 

N. R. Markham and M. Zuker, Nucleic Acids Res. 33, 

W577 (2005). 

B.M. Mladek et al, Multi-scale, ab-initio coarse-graining 
of DNA-functionahzed colloids, in preparation. 
Of course, the discrepancies between many-body inter- 
actions and two-body interactions will be configuration- 
and neighbour-dependent. While the chosen arrangement 
per se is not found in typical crystalls, we chose it for 
convenient comparison with the data of $2- Mechanical 
stability of crystals is mediated by DNA bonds, thus the 
chosen configuration is representative of crystalline ar- 
rangements in that several 'like' next neighbours compete 
for strands of a given 'unlike' colloid. 

D. Gottwald et al, J. Chem. Phys. 122, 204503 (2005). 
J. Fornleitner et al, Soft Matter 4, 480 (2008). 

F. J. Martinez- Veracoechea et al, Phys. Rev. Lett 107, 
045902 (2011). 



